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Abstract 

We study the turbulent velocity dispersion spectra of the coexistent HCN and HCO"'" 
molecular species as a function of length scale in the M17 star-forming molecular cloud. 
We show that the observed downward shift of the ion's spectrum relative to that of 
the neutral is readily explained by the existence of an ambipolar diffusion range within 
which ion and neutral turbulent energies dissipate differently. We use these observations 
to evaluate this decoupling scale and show how to estimate the strength of the plane-of- 
the-sky component of the embedded magnetic field in a completely novel way. 

Subject headings: ISM: clouds — ISM: individual (M17) — ISM: magnetic fields — 
Physical Data and Processes: turbulence 

1. Introduction 

Star formation is one of the very basic topics of astronomy. It was proposed as early as in 
the 18th century, by Immanuel Kant and Simon de Laplace independently, that stars originate 
from clouds of matter contracting under their own gravitation. This point of view is supported by 
modern observations, which show a high correlation between the positions of newly formed stars 
and those of molecular clouds. Molecular clouds are made up primarily of molecular hydrogen with 
a mean number density of approximately 10^ cm~^. Some of the clouds are as big as 10 - 100 pc 
with masses reaching 10^ - 10^ M©, while the typical temperature of an individual cloud is only 
10 - 50 K. Since the Jeans mass, which is the minimum cloud mass needed to overcome thermal 
motion and bring gravitational collapse, is approximately 20 M©, it follows that molecular clouds 
are highly gravitationally bound. On the other hand, the star formation efficiency, defined as the 
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ratio of stellar mass to host cloud mass, is observed to be only a few percents |Myers et alJll986l ). 
Moreover, the free fall time of a molecular cloud (on the order of 1 Myr) is usually shorter than 
its age (approximately 2-10 Myr). These observations imply the existence of mechanisms that 
prevent the rapid collapse of molecular clouds and explain the low efficiency with which stars are 
formed. 

Different mechanisms have been invoked to regulate star formation. If the potential for mag- 
neti c fields embedded in molecular clouds to bring support aga inst self-gravity was recognized early 
on (IShu. Adams. &: Lizand 119871 : iMouschovias &: Ciolekl Il999l ) , the role playe d by interstellar tur- 
bulen ce in the star-forming process was also put forth and is actively studied (jMac Low fc Klessen 
2004 ). The fact that such fundamental questions are still debated owes largely to the lack of observa- 
tional constraints on existing star formation theories. For example, our knowledge o f magnetic field 



stren gths in molecular clouds is mainly based on observa tions of the Zeeman effect (ICrutcher et al 



1993) and the so-called Chandrasekhar- Fermi method ( Chandrasekhar fc Fermi 19531 ). but both 
types of measurements are difficult to carry out in most parts of a cloud. Similarly, turbulent en- 
ergy needs to dissipate at small scales to allow for cloud material to collapse and form stars. But 
the details pertaining to the mechanisms at the origin of turbulence dissipation, whether they are 
related to viscosities or ambipolar diffusion, are still uncertain. 

In this paper, we show how the difference between the turbulent velocity dispersion spectra of 
coexistent ions and neutrals can shed light on the existence of an ambipolar diffusion range, and be 
used to e stimate magnetic field strengths in molecular clouds. To do so, we first use the method put 
forth by lOstriker. Stone. &: Gammid (|200ll ) to study simulations of turbulent molecular clouds to 
measure the spectral index of the Kolmogorov-type turbulence observed in our velocity dispersion 
spectra. We then argue that the observed difference between the ion and neutral spectra is due 
to differences in their respective energy dissipation ranges. Based on the measured spectral index, 
and the difference in the turbulence spectra, we describe how the strength of magnetic fields can be 
estimated. 



2. Turbulence in Molecular Clouds 



The observations of molecular line profiles in most parts of giant molecular clouds (GMCs) 
reveal much wider line widths (on the order of 1 km/s and more) than what is expected from 
thermal motions (a fraction of 1 km/s). While there is no evidence that this motion is due to 



systematic mass motion, e.g., gravitational collapse. IZuckerman &: Evand ([l97J) suggested that the 
large line widths stem from lo calized supersonic turbulence. This idea is strongly supported by 
the discovery by iLarsonl (Il98ll ) that the line width follows a Kolmogorov-type spectral law, i.e., 
a oc L^-^^^ where a stands for the velocity dispersion meas ured from line profiles and L is the linear 
size of the observed clouds or cloud cores (see the Fig. 1 of lLarsonI |l98ll )). Larson's law was derived 
from the statistics of many clouds/cores of different sizes, but the characteristics of the turbulence 
velocity spectrum as a function of the length scale vary from cloud to cloud. This can be seen from 
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Larson's work, where, for example, data pertaining to the M17 molecular cloud and cores embedded 
within it imply a smaller spectral index than that (i.e., 0.38) of the overall law mentioned above. 



Several methods have been proposed to study t urbulence at sub-cloud /core scales (|Ossenkopf &: Mac Low 



2002 ). We apply here the technique introduced by lOstriker. Stone, fc Gammid (|20Qll ) to our data 
obtained for M17. Our data consist of the HON ( J = 4 ^ 3) and HCO+ ( J = 4 ^ 3) rotational 
transitions observed in an approximately 2 arcmin. x 2 arcmin.-wide region located near the center 
of Ml 7 with the Caltech Subm illimeter Observatory (CSO). These data were previously discussed 
in details by lHoude et al.l (120021 ). Figure [1] shows a log-log plot of the HCN velocity dispersion spec- 
trum as a function of the length scale. The angular resolution of the CSO's 10.4-meter telescope 
is approximately 20 arcseconds at the frequencies of interest (i.e., approximately 355 GHz and 357 
GHz, respectively, for HCN and HCO"'" at this transition), and corresponds to the smallest scale on 
the plot. The data at larger scales were obtained by combining the spectra from adjacent regions to 
effectively simulate beam sizes of double, triple, and quadruple the original size. In Larson's work 
it was adequate to approximate the clouds or cores as being spherical in shape, and scales along the 
line of sight could, therefore, be reasonably estimated from corresponding spatial extents on the sky. 
On the other hand, at sub-cloud/core scales the sky-projected and line-of-sight scales are no longer 
directly correlated. Consequently, the scatter at eac h scale in our data of Figure [J is mainly due 
to the range of line-of-sight scales. But, interestingly. lOstriker. Stone. Sz Gammid ([200l|) showed in 
their simulations of turbulent molecular clouds that the lower envelope of such velocity dispersion 
distributions traces the actual turbulent velocity spectrum very well. 

Two points are worth noting. First, the lower envelope of our HCN data shows a spectral 
index o f 0.18 (as measured from the logarithmic slope), which is smaller than that observed by 
Larson (1981) for his ensemble of molecular clouds. Second, in contrast to the simulations of 



Ostriker. Stone, fc: Gammid |200ll ). this spectral index is constant acr oss the dispersion spectrum. 
The c hange in spectral index at smaller scales seen in the results of lOstriker. Stone. &: Gammie 
|200ll ) (see their Figure 10(a)) is explained by the presence of a turbulent energy dissipation scale. 
Accordingly, we interpret the constant spectral index of the HCN lower envelope in our data as an 
indication that we did not reach the dissipation scale in Ml 7 at the spatial resolution with which 
our data were obtained. That is, our data probe the inertial range. 



3. The Ambipolar DifTusion Scale 

It is usually believed that magnetic fields in the interstellar medium are highly coupled to the 
gas; this is the so-called "flux freezing" phenomenon. This is one of the most important consequences 
of the induction equation 

— = Vx (vxB)+77V2b, (1) 
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where B is the magnetic field, v the velocity of the (ionized) gas, and r] the magnetic diffusivity 
(assumed to be constant here). The ratio of the convection term to the diffusion term is defined as 
the magnetic Reynolds number 



|V X (v X B) 
|?7V2B| 



VL 



(2) 



where L and V are characteristic spatial and velocity scales, respectively. Under condi tions when 
Rm ^ 1, the flux freezing approximation is warranted through Kelvin's vorticity theorem (jChoudhuri 
19981 ). as the convection term dominates over the diffusion term in equation ([T]). 



When dealing with weakly ionized plasmas, however, the diffusivity is not limited to that 
conveyed through rj. The corresponding analysis reveals a more general form for Ohm's law (jParker 
2003) 



E = l 

c 



„ 47r / . . B ^. 

-Vn X B + — 77J + aj X — + 
c V ^ 



(3) 



where c is the speed of light, Vn is the velocity of the neutral component of the gas, j and are 
the current density and its component perpendicular to B, while a and f3 are, respectively, t he Hall 
coefficient and the effective magnetic diffusivity (also known as the Pedersen coefficient (jParker 
20071 )). It can be shown that for typical conditions encountered in GMCs P ^ a ^ tj, and the 
above relation for the electric field thus simplifies to (lBiskampll2003l ) 



E 



-- I Vn X B 



c 



(4) 



The effective magnetic diffusivity is given by 



47rni/xt'i ' 



(5) 



where n\ is the ion density, while v\ and // are the collision rate of an ion with the neutrals and the 
mean reduced mass characterizing such collisions, respectively. Although the insertion of equation 
(111) into Faraday's l aw of induction does not viel d a relation analogous to equation ([T]) for the 
induction equation (IBrandenburg h Zweibell 1 19941 ) , it is nonetheless clear that the flux freezing 
approximation will cease to be adequate when the effective magnetic diffusivity term dominates 
over the convective term in equation ([4]). Consequently, we can define a new effective magnetic 
Reynolds number with 



Rn 



C |Vn X Bj 
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where Ampere's law was substituted for the current density. Based on equation ([6]), it is clear that 
although Rra 3> 1 for scales usually probed with single-dish observations (i.e., on the order of sizes 
subtended by the telescope beam), this relation does not necessarily hold true for all the turbulent 
eddies existing in any given region of a molecular cloud. This is because situations involving small 
eddies (in the neutral gas) of size L evolving in a region embedded with a (strong) magnetic field 
will lead to lower effective magnetic Reynolds numbers. In fact, it is expected from the loss in the 
relative importance of the convective term in equation (|4| that the magnetic field, and therefore 
the ions, will decouple more easily from neutral eddies at such scales, when the effective Reynolds 
numb er is on the order of unity or less. This decoupling, or ambipolar diffusion (IMestel &: Spitzer 
19561 ). will be responsible for friction forces between the ion a nd neutral components of t he gas, 
and will thus bring about the dissipation of turbulent energy (jZweibel &: Josafatssonlll983l ). It is 
also reasonable to expect that the existence of ambipolar diffusion will be the cause for possible 
differences between the velocity dispersion spectra of the ions and the neutrals. 

We also note that for typical conditions in the interstellar medium (e.g., B ~ 10 ^G, ^ 10^ 
cm~^, an ionized fractio n of approximately 10~^, and with Ui ~ 1.5 x 10~^nn s~^ for ion/neutral 



collisions (INakand 11984 )) equation jS]) yields an effective magnetic diffusivity /3 of the order of 
10^^ cm^s~^, which is much larger than the hydrodynamic viscosity (which is of the order of 10^^ 
cm^s~^). The situation also holds for denser parts of molecular clouds such as the ones probed 
with our data, where both density and magnetic field strength can increase by as much as two 
or three orders of magnitude. Given this fact, will ambipolar diffusion set an energy dissipation 
scale in the turbulent energy spectrum that is larger than the scale due to viscosity? This is 
still an open question, both on observational and theoretical/simulation grounds. For supersonic 
magnetohvdrodvna mic turbulence, the fast waves and Alfven waves are strongly damped at smaller 



scales when < 1 ( McKee et al 



1993 



Balsaralll996l ). It is, therefore, reasonable to expect steeper 



energy spectra at ambipolar diffusion scales. But since ions and neutrals are decoupled from each 
other, their respective energy (and velocity dispersion) spectra must have different slopes at these 
scales. Which one should be steeper? Here we show that observing the turbulent velocity spectra of 
coexistent ion and neutral molecular species, even at spatial resolution much lower than ambipolar 
diffusion scales, can shade some light on these questions. 



4. Comparison of HCN and HCO+ in M17 

Figure [2] shows in part the same data as in Figure [H with the exceptions that the vertical 
axis is now for the square of the velocity dispersion, we used linear scales for both the horizontal 
and vertical scales, and more importantly data for the other molecular species, HCO"'", have been 
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added. The use of linear scales for the axes in this figure was chosen to better show the difference 
between the two data sets as a function of the length scale, as this will become important for the 
discussion below. As can be seen, the two molecular species show very similar spectral indices for 
their respective velocity dispersion spectra (i.e., 0.18, as calculated with the lower envelopes), which 
is expected at this scale because ions are trapped by magnetic fields, which are in turn well-coupled 
with neutrals due to flux freezing. But note that the lower envelope of HCO"'" is downshifted from 
that of HCN. We now turn our attention to identifying the cause for this nearly constant shift of 
the ion velocity spectrum over all the scales shown in the part of the inertial range probed by our 
observations. 

HCN and HCO"'" hav e been showii to be c oexist e nt in molecular clouds and to have system- 
atic line width differences (IHoude et al.ll2000; 



12002 



Lai. Velusamy. fc Langer 



2003 



Houdeetal 



2004 ). Due in part to its smaller critical density HCO"*" is generally somewhat more optically 
thick and will accor dinglv usually exhibit slightly more extended spatial distributions than HCN 
iHoude et al.ll2000al |hlv Both of these properties would suggest that HCO"*" should have wider line 
widths than HCN, yet the fact is t hat its spectra l line p rofiles are consistently narrower in turbulent 
molecular clouds (such as M17). iHoude et al.l (|2000al ) put forth a model, where magnetic fields 
affect the dynamics of molecular ions through the cyclotron interaction, that can account for the 
narrowin g of their spectral l ine profiles, which was in turn interpreted as a signature of ambipolar 
diffusion iHoude et all hooi ) . 



Here we further argue that the aforementioned observed difference between the two turbulent 
velocity spectra (as a function of length scale) is also the result of ambipolar diffusion. Even though 
ambipolar diffusion happens at much smaller scales than that at which our data was obtained, it 
is still discernible in our observations because the velocity dispersions we measure stem from the 
contribution of all turbulent eddies with scales approximately equal to, or smaller than, the different 
beam sizes. To explain the systematically downshifted velocity dispersion spectrum of the ions, it 
is necessary that 



1. Ambipolar diffusion does not set the cut-off range for both the ion and neutral energy spectra. 

2. The ion spectral slope must be steeper than that of the neutrals' at ambipolar diffusion scales, 
as illustrated in Figure [3] (more below). 



As mentioned before, a steep drop in the ion velocity dispersion spectrum at ambipolar diffusion 
scales can be expected due to the fact that most of families of magnetohydrodynamic waves are 
strongly damped when Rm < 1- Neutrals, while still partaking in the turbulent energy cascade 
of the inertial range, should also exhibit a steeper turbulent spectrum than in their inertial range 
because of friction forces between ions and neutrals. But since the low ionization level will limit the 
amount of friction, neutrals will still have their energy cut-off set by hydrodynamic viscosity. The 
square of the observed velocity dispersion at a particular scale is proportional to the integral of the 
energy spectrum below that scale. Therefore, at all scales larger than those covered by the ambipolar 
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diffusion range the difference between the squares of the neutral and ion velocity dispersions will 
be proportional to the difference of the corresponding integrals between the scale where ambipolar 
diffusion sets in and the hydrodynamic dissipation scale (see below). This is conveyed through the 
shaded region in Figure [3l 

5. Estimation of the Magnetic Field Strength 

The difference between the ion and neutral velocity dispersion spectra can be used to estimate 
the strength of the magnetic field. To see how this can be achieved, let us consider the wave number 
k' = 2t\ jV corresponding to turbulent eddies of characteristic size V such that = V^L' / j3 = 1, 
with the characteristic velocity of these (neutral) eddies. As discussed earlier, this is the scale 
at which ambipolar diffusion sets in. For a wave number K < k' (i.e., of spatial scale larger than 
L'), the difference between the squares of the neutral and ion velocity dispersions can be written as 

/■oo 

= / [Fn (k) - Fi (k)] dk, (7) 
Jk 

where F^{k) and Fi(k) are (twice) the turbulent energy spectra density (i.e., per unit mass and 
wave number) for the neutral and ion species, respectively. But because of the good coupling of the 
ions and neutral through flux freezing in the inertial range we expect that 

Fn{k) =Fi{k), fork < k' . (8) 

If we further assume that Fi (k) decreases steeply beyond A;', then for any wave number K < k' 
the difference of the square of the velocity dispersions (eq. [7]) is only a function of k' with 

Aal. ~ Aal = Jj^ (k) dk, for K < k' . (9) 

This result is consistent with our observations, since, as was mentioned earlier, our data show 
such a nearly constant difference between the neutral and ion spectra. More precisely, the measured 
differences Aaj^ from the velocity dispersion spectra of Figure [2] are 0.57, 0.65, 0.57, and 1.22 km^/s^ 
with increasing beam sizes. The discrepancy at the largest spatial scale is easily accounted for by 
the small amount of statistics available, which prevented us from finding the true lower envelope 
for both velocity dispersion spectra. 

An important consequence that follows from this analysis is that the determination of the 
decoupling scale L' and the neutral velocity dispersion at L' would allow for the evaluation of 
the magnetic field strength through equation ^ by setting = 1. More precisely, we find that 



52 = ATTUiflUiV^L' , 



(10) 
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which can be conveniently transformed using typical parameters for GMCs. For example, if we use 
mean neutral and ion masses of, respectively, 2.3 and 29 times that of the hydrogen atom, then 

p ( L' y/y K V^V \( xi y/^ ^ .... 

V0.5mpcy Vlkms-V V 106 cm-3 / V lO-^/ ' ^ ' 

where Xi is the ionization fraction. 

It is important to note that since observations can only reveal the velocity dispersion along 
the line-of-sight and that the effective magnetic diffusivity term in equation ([IJ is in a direction 
perpendicular to the magnetic field, the strength of the field estimated with equations (fTO]l or ([TT]) 
is only valid for its plane-of-the-sky component. 

We now proceed and evaluate L' and using our data. First, we note that the velocity 
dispersion due to an eddy of scale size k = 2tt/L \s given by 



(k) = (k) Ak. (12) 
The energy density (k) can readily be evaluated from our data, since we know that 



/•oo 

(^n (^) = I Fn 

Jk 




where (k) is the neutral velocity dispersion, while b and n are the different fit parameters 
to the neutral data (see Figure [2]). Using the derivative relative to k on this equation and the 
spectral spread corresponding to a Gaussian beam of full-width-half-magnitude L to estimate Ak 
(= Y^Bln (2)/L ), we get from equation (fT2]l that 

(L) ct 0.37bnL'' . (14) 

Moreover, if, as assumed, there is a significant drop of the ion spectrum in its dissipation range 
(as shown in Figure [11), it then follows that the ion velocity dispersion (Ji{k') observed at L' (or 
k' = 2t: /L') will be approximately equal to F^. We, therefore, write 

+ (15) 

where it should be clear that —a corresponds to A(t|,, as defined in equation ([9]) and measured from 
our data (see Figure [21). Equating equations (fT4ll (at L') and (fTSl) gives 
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For M17, the fits to our data yield a ~ — 0.59km^s~^, b ~ 1.17km^s~^arcsecond~"^, and 
n ~ 0.36. Equation (fT6]) then gives L' ~ 0.22" or 1.8 mpc at the distance of M17 (i.e., approximately 
1.7 kpc), and ~ 0.30 km/s. 

A precise determination of the magnetic field strength in Ml 7 through equation (fTTI) would 
require an equally precise knowledge of the gas density rin and the ionization fraction Xi, which we 
do not possess at this time. We can, however, use reasonable values for these parameters and provide 
an estimate of the magnetic field strength to allow for comparison with previous measurements 
obtained with other techniques. Therefore, choosing Un = lO^cm"^, which should be appropriate 
for HCN and HCO"'" at the J = 4 ^ 3 transition, and Xi = 10~^ for the ionization fraction we 
calculate a plane-of-the sky magnetic field strength of B ^ 1.0 mG from equation (fTTI) • Although 



this value compares well with earlier HI and OH Zeeman measurements (jBrogan fc Trolandl 120011 ) 
that yielded line-of-sight magnetic field strengths as high as 0.75 mG in this molecular cloud, we 
should keep in mind that our estimate is probably uncertain to within a factor of a few or even as 
much as an order of magnitude. We also note that the details of our calculations rest in part on 
the assumption that the ion velocity dispersion spectrum drops steeply in the ambipolar diffusion 
range. In the event that this assumption is not realized, the characteristic velocity V^, and the 
magnetic field strength, would be overestimated through equation p^ . 



It is important to note that current (e.g., SMA, CARMA) or future (e.g., ALMA) interfer- 
ometers could provide high enough spatial resolutions to allow for the direct determination of the 
ambipolar diffusion scale L' and the characteristic velocity of an eddy at that scale directly from 
the velocity dispersion spectra. The determination of L' would be straightforward from the velocity 
dispersion spectrum of the ion, as it would be identified with the scale at which a significant drop 
in velocity occurs (i.e., the ion cut-off point) or where the neutral and ion spectra depart from one 
another (see Figure [3|). 



6. Summary 

Our analysis of turbulent velocity dispersion spectra from coexistent ion and neutral molecular 
species can be summarized as follows: 



1. We showed th at by combining observation a l spec tral line profile data to different scales, in 
the manner of lOstriker. Stone, fc Gammid (j200l|), it is possible to probe turbulent velocity 
dispersion spectra at sub-cloud/core scales. 



2. For coexistent ions and neutrals in M17, the turbulent velocity dispersion spectra share the 
same slope in the inertial range, as expected from the flux freezing approximation, but the 
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ion spectrum is downshifted nearly uniformly at all scales relative to that of the neutrals, as 
shown in Figure [2l This can be explained by the occurrence of ambipolar diffusion at smaller 
scales where i?m < 1. The downshifted velocity spectrum of the ions implies that they have 
a steeper turbulent energy spectrum than the neutrals at ambipolar diffusion scales. It is, 
therefore, apparent that ambipolar diffusion does not set the energy cut-off range for both 
ions and neutrals. 

3. Since ambipolar diffusion happens in the presence of magnetic fields, the difference between 
the velocity dispersion spectra of ions and neutrals can be used to estimate magnetic field 
strengths. For M17, we estimate the ambipolar diffusion scale to approximately 1.8 mpc, 
and the plane-of-the-sky magnetic field strength to approximately 1 mG. Because we do not 
precisely know the gas density and the ionization fraction in M17 this estimate of the magnetic 
field strength is probably uncertain to within a factor of a few or even as much as an order of 
magnitude. 



We note that using optically thin molecular species, such as the H^^CN and H^^CO"^ isotopologues, 
co uld be advantageous for t his type of study; their sy s temat ic line with difference is also reported 
by iHoude et al.l (l2000bl ) and lLai. Velusamy. fc Langen (|2003l ). Their advantage would reside in the 
fact that their corresponding line profiles would not be saturated, unlike the species used in this 
work in some parts of the molecular cloud under consideration, and the line widths measured would 
reflect more faithfully the true velocity dispersions of the ion and neutral gas components along 
the line-of-sight. Moreover, H^^CN and H^^CO^ would allow us to probe magnetic fields at gr eater 
depth in star-forming regions, where they are likely to be stronger |Crutcheiill999l : lBasull2000l ). 
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Fig. 1. — The HCN velocity dispersions a obtained from spectral line width measurements at sub- 
cloud scales L (using logarithmic scales). The lower envelope can be fitted to a Kolmogorov-type 
spectrum with a spectral index of approximately 0.18, as measured from the logarithmic slope. The 
observed data are at the smallest scale shown (i.e., 20 arcseconds), while adjacent positions were 
combined to get velocity dispersions at larger scales. 
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Fig. 2. — Similar to Figure [H with the exceptions that the vertical axis is now for the square of 
the velocity dispersion, we used linear scales for both the horizontal and vertical scales, and the 
data for the other molecular species, HCO"'", have been added (in red, HCN is in black). The HCN 
data are fitted for a Kolmogorov-type law = hL^ to their lower envelope. The relative downshift 
of the HCO"'" velocity dispersion spectrum is made explicit with the corresponding a? = a + bL'^ 
fit. The different parameters are a ~ — 0.59km^s""^, h ~ 1.17km^s~^arcsecond and n ~ 0.36. 
The measured differences between the two velocity dispersion spectra are 0.57, 0.65, 0.57, and 1.22 
km~^s~^ with increasing beam sizes. The lack of agreement at the largest spatial scale is easily 
accounted for by the small amount of statistics available, which prevented us from finding the true 
lower envelope at that scale for both spectra. The uncertainties on the measured velocity dispersions 
for the lower envelopes are negligible (< 0.01 km/s) and are thus not shown. 
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Fig. 3. — An illustration of our explanation for the difference between the ion and neutral velocity 
spectra of Figure [2] is shown. This is a log-log plot of the turbulent energy versus wave number 
(inverse of the spatial scale) in arbitrary units. Ambipolar diffusion happens at wave numbers larger 
than k' , which is in turn smaller than k" , where the hydrodynamic viscosity sets in for the neutral 
spectrum (k) (Fi (k) is for the ion spectrum). The energy spectrum in the inertial range, F {k), 
is common to both species, while ions and neutrals have different spectra at ambipolar diffusion 
scales. The square of the velocity measured at a particular wave number K is proportional to the 
integral of the energy spectrum over all wave numbers greater than K. The observed difference 
between the two velocity dispersion spectra is proportional to the shaded area. 



